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1 Introduction 



Gauge theories in higher dimensions are a promising candidate beyond the Standard 
ModeL Such theories turn out to possess unexpectedly rich properties that shed new 
hght and give a deep understanding on high energy physics. In fact, it has been shown 
that new mechanisms of gauge symmetry breaking [1, 2, 3, 4, 5], spontaneous supersym- 
metry breaking [6], and breaking of translational invariance [7, 8] can occur, and that 
various phase structures arise in field theoretical models on certain topological manifolds 
[9, 10, 11]. Furthermore, new diverse scenarios of solving the hierarchy problem have been 
proposed in [12, 13, 14, 15]. 

Since the origin of gauge symmetry breaking is still an unsolved problem, it is worth 
pursuing an alternative mechanism to mimic the Standard Model Higgs. In this paper, 
we focus on the gaugc-Higgs unification that the "Higgs" field arises from an extra di- 
mensional component of a higher dimensional gauge field [1, 2, 4]. Higher dimensional 
gauge invariance forbids any Higgs mass term at tree level and the Higgs potential can be 
generated through quantum corrections. In the gauge-Higgs unification, the Higgs field 
corresponds to the Wilson line phase and the effective potential never suffers from ultravi- 
olet (UV) divergences due to the nonlocal property of the Wilson line phase. As a result, 
the vacuum expectation value and the mass of the Higgs field derived from the effective 
potential are finite and calculable. This is a very attractive feature of the gauge-Higgs 
unification. A fiat compactification on a circle S^, however, leads to the light Higgs mass 
problem in normal settings of gauge-Higgs unification models [16]. A compactification 
on a warped extra dimension may solve the problem [17]. In this paper, we take an 
alternative approach to solve it. 

In the standard gauge-Higgs unification, the kinetic terms of bosonic (fermionic) fields 
in the action are taken to be the second (first) order derivatives for extra dimensions as 
well as the Minkowski space-time. Lorentz invariance does not, however, require the order 
of the derivatives to be the same for the extra dimensions and the Minkowski space-time 
because it is explicitly broken by the compactification. We may thus have an opportunity 
to introduce higher derivative terms for the direction of extra dimensions. 

Recently, Horava [18] has proposed an interesting idea to make gravity theory power- 
counting renormalizable in 4-dimensions. His idea is to treat space and time non-relativistically 
and the anisotropy between them is given by the introduction of higher spatial deriva- 
tive terms characterized by the dynamical critical exponent z. A number of studies have 
been made on the Hofava-Lifshitz gravity.^ Field Theories with higher spatial derivative 

''studies at an early stage have been given in Ref.[19] for cosmological applications, Ref.[20] for black 
hole physics and Ref.[21] for theoretical aspects. 
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terms have also been investigated on various subjects of renormalization [22], non-gauge 
models [23], QED [24, 25], Yang-Mills theory [28], SUSY [29, 30] and the Standard Model 
extension [31]. All Hofava-Lifshitz type models lose Lorentz invariance at high energies 
but it is expected that it would emerge at low energies as an accidental symmetry [18]. 
The recovery of Lorentz invariance at low energies is, however, a nontrivial problem since 
in a theoretical point of view there is no reason that different particles possess the same 
limiting speed (i.e. the common light velocity c) in the absence of Lorentz invariance.^ 
Problems of Hofava-Lifshitz type theory have been reviewed in Ref.[32]. 

In this paper, we introduce higher derivative terms only for the direction of extra 
dimensions and keep Lorentz invariance intact for the Minkowski space-time. Thus, our 
models have the anisotropy between the Minkowski space-time and the extra dimensions 
but not between space and time. This situation will be rather close to the original 
idea proposed by Lifshitz [33]. Higher order derivative terms would become important 
if coefficients of lower order derivative terms happen to vanish. This is indeed the case 
when the system lies at a Lifshitz point [34] . We thus consider the gauge-Higgs unification 
in a higher dimensional gauge theory at a Lifshitz point, though we will not investigate 
whether our models would lie at a Lifshitz point, but simply assume it in this paper. 

As noted before, the effective potential for the Higgs field, which is a zero mode of 
an extra dimensional component of a higher dimensional gauge field, is finite and free 
from UV divergence. This UV insensitivity does not, however, imply that higher order 
derivative terms are irrelevant to the effective potential. It turns out that they bring 
about quantitative changes of the effective potential and play an important role to solve 
the hght Higgs problem. To show this is one of the main purposes of this paper. We will 
also develop some mathematical tools to compute one-loop effective potentials. This is 
another purpose of this paper. 

This paper is organized as follows: In the next section, we explain our setup of the 
gauge-Higgs unification in a framework of Lifshitz type gauge theory. In section 3, we eval- 
uate the one- loop effective potential for the zero mode of the extra dimensional component 
of the gauge field and discuss how higher derivative terms affect the effective potential. 
In section 4, we present a five dimensional SU{2) model to demonstrate properties found 
in section 3, explicitly. In section 5, we extend the results in section 3 to the Lifshitz type 
gauge theory on R^~'^ x x S^. Section 6 is devoted to conclusions and discussions. 
Technical details will be found in appendices. 

resolution to this problem has been proposed in Ref.[30]. 
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2 Lifshitz Type Gauge Theory on R^-^ X 



In this section, we present our setup of the gauge- Higgs unification in a framework of 
Lifshitz type gauge theory. To this end, we consider a D-dimensional SU{N) gauge 
invariant theory compactified on a circle.^ 



where x^^ {fj, — 0,1, D — 2) denotes the {D — l)-dimensional Minkowski space-time 
coordinate and y is the coordinate of the extra dimension on the circle of the cir- 
cumference L. The covariant derivatives and Dy for the field ^ belonging to the 
representation TZ of SU (N) are given by 



where T^'^^ (c = 1, 2, ■ ■ ■ , N"^ — 1) is a generator of SU{N) in the representation TZ. 

We should here make several comments on the higher derivative terms in the action 
(1). We have introduced the higher derivative terms only for the direction of the extra 
dimension and assumed that there are no higher derivative terms such as {d^d'^)"^{dy)"'. 
This implies that there is no possibility to add any higher order derivative terms for the 
gauge field in a gauge invariant way because Fyy = and the term {F^yF^^Y for z > 1 
would produce higher derivatives with respect to the extra dimensional coordinate and 
the Minkowski space-time ones as well. We could add lower order derivatives of the extra 
dimensional one such as (i^yDy)"' {n = 1,2, - ■ ■ ,z — 1) for the fcrmion and {DyY^ for 
the scalar, but we have omitted those terms because we would like to clarify the effects 
of the higher derivative terms. It could be justified if the system sits on a Lifshitz point, 
where the lower derivative terms become irrelevant. We will not, however, pursue such 
a possibility. We simply keep only the highest order derivative terms labeled by the 
dynamical critical exponent z and investigate the theory at a Lifshitz point. It should 
be noticed that a length parameter / has been introduced to adjust the dimension of the 
higher derivative terms. The ratio l/L turns out to be very important in determining the 
vacuum expectation value and the mass of the Higgs field. 

Since the extra dimension is compactified on the circle, we have to specify boundary 
conditions on the fields. We here take the periodic boundary condition for all the fields, 

^ Since the gauge fixing and the ghost terms are irrelevant in our discussions, we have omitted them 



s 




(1) 



Dym^{dy + igAlTlT^))m, 



(2) 
(3) 



in Eq.(l). 
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A^{x,L) = A^{x,0), 
Ay{x,L) = Ay{x,0) , 

$(x,L) = $(x,0). (4) 

The extension to other boundary conditions will be straightforward. 

In the gauge- Higgs unification, the zero mode of Ay plays a role of the Higgs field and 
the vacuum expectation value {Ay) can be determined dynamically as the minimum of 
the effective potential, which is induced through radiative corrections. The purpose of 
the next section is to compute the one-loop effective potential for the zero mode of Ay 
from the action (1) with the boundary conditions (4). 

3 One-Loop Effective Potential 

In this section, we evaluate the one-loop effective potential for the dynamical variable 
(Ay) from the action (1) of the 5d Lifshitz type gauge theory on R^"^ x S^. 

Kff («) = Scalar Ms) + V^fermion Mp) + V^gauge ' («) > (5) 

where the subscript z of indicates the Lifshitz type higher derivatives of the dynam- 
ical critical exponent z for the direction. V^faiar ^^^^ ^ermioiT^'' ^auge ^'^^ denote 
the contributions from scalar, fermion and gauge loop diagrams, respectively and the 
dimensionless variable a is defined by 

a^g{Ay)^. (6) 

In the following, we first evaluate the effective potential for the massless matter and then 
for the massive one. 

3.1 Massless Matter 

In this subsection, we consider the massless scalar and fermion fields with Ms — Mp — 0. 
Let us start with Vg^aiar K'^''> ^s — 0) which comes from the scalar loop diagram 

Kcaiar ' {a; Ms = 0) = J - 2^ tr(^) Ini^PE+l"- \ ^ ) ) 

m=—oo '- 

(7) 



^Although the periodicity of the fields on a circle may require the boundary conditions (4) up to 
gauge transformations, we here restrict a class of the gauge parameter K{x,y) to the periodic function 
A(x,L) = A(a;,0), so that all the fields are assumed to be periodic with respect to the coordinate y on 
S\ 
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where denotes the (D — 1) -dimensional Euchdean momentum. When the Lifshitz scalar 
$ which propagates the loop belongs to the representation TZ of the gauge group, the trace 
in Eq.(7) should be taken over the gauge indices with respect to the representation TZ and 
a should be expressed os a — a'^T^^y We note that Vgfaiar unrenormahzed 
quantity, so that it should be renormalized to obtain a finite expression. 

For z — 1, which corresponds to the normal kinetic term, the techniques have already 
been well estabhshed to obtain the finite expression from Eq.(7). Some of the techniques 
for z = 1, however, turn out not to apply for the case of 2; > 1. Thus, we need to develop 
mathematical tools to compute the one-loop effective potential for 2; > 1. This is one of 
the purposes of the present paper. 



Let us now evaluate VI 



scalar 



. We first rewrite Eq.(7) into the form 



scalar 



{a; M5 = 0) = - 



m=— 00 



P+l 

dtt~— 



X tr 



(7^) 



exp{-if(^-i)( 



/27r(m + q;)\ 2 



where we have used the standard formulas 

/ 1 

In ^ = — - 



ds\r{s)J, 



dtf-'e-^' 



d / r 



ds vr(s) . 

dp 



s-)-+0 

A-Kt 



-tp 



For 2; = 1, we can use the Poisson summation formula 

'27r(m + a)\2^ L 



exp{-t(- 



} 



(8) 

(9) 
(10) 

(11) 



(12) 



and then arrive at the well known renormalized expression^ 

T{D/2) 1 



KfJa:''^^'^H«;^5 = o) = -2: 



ID 



n=l 



n) 



cos{27ma) 



n 



D 



For 2; > 1, we may use the original Poisson summation formula, i.e. 

, , J —00 



—2-Kinu 



(13) 



(14) 



^The UV divergent term of n = has been removed from the summation in Eq.(13) to renormahze 
the one-loop effective potential. 
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mplaceofEq.(12). With the help of the above formula and the relation r(s) = dtt^ *, 
we can have a renormalized expression 



1 \ ^ 

i \ 2 



X 



^'tr(7^)[exp{27rma}] / du\u\'^^-^^e-^'''"'' , (15) 

^ J —OO 



where the prime of the summation denotes that the n = contribution, which corresponds 
to the UV divergent part, should be removed from the summation to get a finite result of 
the effective potential. Furthermore, using the formula of the Fourier transform 



/_ 



2 / |27rn|i+^ 



(16) 



for i/ 7^ • • • , —3, —1, 0, 2, 4, • • • , we obtain 



Kfala"''''(«;M5 = 0) 



D + l 
2 



X 



1 

4:71 

l{z-l){D-l) ^ 
n=l 



V{-{D - l)/2)r(l + ;z(D - l))sin 
cos(27rna) 



^^7^^CD-l)^j 



n 



D+(^-l)(£)-l) 



(17) 



It should be noticed that the expression (17) would be ill defined for D — odd or 2; = even 
because the formula (16) cannot apply for those cases and also the factor V{—{D — l)/2) 
diverges for D — odd. To avoid these problems, we regard the space-time dimension D as 
a real (or complex) number and define the expression (17) by the analytic continuation 
of D. It then turns out that the final result is finite and well defined for any positive 
integers D and z: 



1 . ; . (^-l)(D-l) 



cos(27rnQ;) 



n=l 



n 



D+(z-l)(D-l) 



(18) 



where 



sin(7r2:(s-l)/2) 

hn z = hm — ^ , / 

' s^D sin(7r(s- l)/2) 



z for D = odd, z = odd, 

^(_l)(i^-i)/2 for D = odd, z = even, 

(_l)(^-i)/2 for ^ even, z = odd, 

ior D — even, z — even. 



(19) 



Here, we have used the formulas of the gamma function 



Sm TTS 



i2s 



r(2.) = 2^r(.)r(s + i/2). 



(20) 
(21) 
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For z — 1, Eq.(18) is found to exactly reproduce the expression (13). This will give a 
consistency check of our result (18). 

To confirm the validity of the result (18) furthermore, we would hke to derive Eq.(18) 
in a different way. To this end, we rewrite Eq. (8) as 



r(-iZ-l) /27r\^(^-') 



-z{D-l) 



V "/ ^ m=— oo ^' 

Then, using the Hurwitz zeta function 

oo ^ 

C(g,a) = ^ for 0<a<l, (23) 



m=0 ^ ' 



we have 



D-l 
2 



D-l 
2 



(47r)^ 

xtr(^){C(-^(L' - 1), a) + C{-z{D - 1), 1 - a)}. (24) 



The formula [36] 

2r(l — s) f . /7rs\ cos(27mQ;) /7rs\ sin(27rnQ;) "i 



for Re s < 0, leads to the same result (18), as it should be. It is interesting to note that 
UV divergence has been removed automatically from the effective potential thanks to the 
zeta function regularization. We will further find other different derivations of Eq.(18) in 
Appendix A. 

We have succeeded to compute the scalar contribution to the one-loop effective po- 
tential for the zero mode of Ay. It is now easy to obtain the fermionic contribution 
^felion = 0). The result is 

' / / \ ^ r cos(27rna) ] , , 



n=l 



The difference between the scalar and fermion contributions appears in the overall factors. 
The factor —2 in V^sf^iar ^ '''^ replaced by -|-2[^/^] in V^f^mion^"^^ . The difference of the overall 
sign comes from that of the spin and statistics for bosons and fermions. The factor 2 in 
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^faiar ^^'^ 2^^^'^^ in V^f^^^^^^'' are just the number of (on-shell) degrees of freedom for 
a complex scalar and a Dirac spinor, respectively.^ 

The gauge field contribution to the one-loop effective potential in our model is nothing 
but that of z = 1, i.e. 

T.J?^-ix5i/ N _ (r, 1 ^ r cos(27rna) i 

^gauge n«) - -(^ - 2)— 2^tr(^dj) -5 • (27) 



n=l 



This result will lead to an interesting observation to determine the vacuum expectation 
value of a. Since V^auge ^ is independent of I, the gauge loop contribution will become 
less important than the scalar and fermion loop ones when the ratio l/L is large. On the 
other hand, it becomes important when l/L is small enough. 

Before closing this subsection, we would like to make comments on the ^-dependence 
of Scalar fermion ■ ^ — ^1 ^hc infinite Summation over n is given by the form 

5^^cos(27rnQ;)/n^. The power of n is just the total space-time dimension D. For general 
z, the power of n is replaced by D + {z — 1){D — 1). This can be understood from an 
anisotropic scaling point of view. It follows from the action (1) that we may introduce 
the anisotropic scahng for the space-time coordinates as [18] 

x" ^ b^x^", (/x = 0,l,2,--- ,D-2), 
y ^ by. (28) 

Then, the effective dimension D of the system may be defined from the scaling of the 
volume element as 

[d^'-^xdy] ^ b^id'^-'xdy] (29) 

with D = z{D — 1) -I- 1 = (2; — 1)(D — 1), which agrees with the power of n in Eq.(18) 
and Eq.(26). 

Another interesting observation is the overall sign of V^^^^y^^ ^ and Vf^^.^^^^ \ For 
z = 1, the overall sign is determined by the spin and statistics but not the space-time 
dimension D. In fact, h^^z = 1 for z = 1 (and any D) but it can change the sign for z > 1 
without changing the spin and statistics. Furthermore, the contribution to the one-loop 
effective potential vanishes for z = even and D = even.^° This is not the case for z = 1. 

The other important observation is that a huge numerical factor can arise for z > 1. 
The numerical values of r(D + (z - 1){D - l))/r{D) for L> = 5 and ^ = 1, 2, 3, 4 are, for 



^For a real scalar or a Weyl/Majorana spinor, the factor should be replaced by 1 or 2[-^/^l ^. 
similar property has been found in the Casimir effect [25]. 
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instance, listed below: 

( 1 for 2; = 1, 

r(F> -i- A( — ^\\ I i7vin3 f^y^ — o 

(30) 



1 for 2; = 1, 

r(5 + 4(^-1)) I 1.7 X 10^ for z = 2, 

r(5) 1 2.9 X 10^ for z = 3, 

8.7 X 10^^ for z = 4. 



Another source of a large/small number will come from the factor (//L)^^-^)^^-^), which 
depends on the ratio of I and L. In the gauge-Higgs unification, the Higgs mass squared 
will be given by 



2 5Veff 
777," ~ 



(31) 



2/ 

Those numerical factors can change the value of the Higgs mass drastically.^^ Therefore, 
the overall factor of V^efr is very sensitive to the magnitude of uih- The Higgs mass 
prediction in Lifshitz type gauge theory is expected to be quite different from ordinary 
gauge-Higgs unification models with z — 1. 

3.2 Massive Matter 



In this subsection, wc evaluate the one-loop effective potential for massive matter. The 
computations of Kjcaiar ''{'^'j^^s) and V^^^j^-^^^^ ^ {a; Mp) will involve a number of tech- 
nical manipulations as compared to the massless case. We here give only the result of 
K^^,:,"'''H^; Ms) ioTD = odd: 

z z{2q+2)-l ^ 

xtr(^) J]Re[(-za;,)^(^^+^) J] ^ Li,(2,+2)-.+i(e-''^+^^-'^)] , (32) 



k\ 

j=l k=0 



where 



ujj = exp(i ^-^^^ ^ tt) , (33) 

l/z 



^'s{x)-J2-. (35) 



n=l 



^^It has also bcicn reported in [26] that the Higgs mass can be heavy due to the effect of the violation 
of the fivc-dinicnsional Lorentz invariance. 
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The details will be found in Appendix B. For D — 5 and z — 1, Eq. (32) reads 

V^2fHa;Ms) = -^tr(^)Re[3Li5(e-^^-+'2-) 

+3LMs Li4(e-^^«+''^") + L'M| Li3(e-^^«+''^'^)] , (36) 

which is consistent with the result in Ref. [37]. It is not difficult to verify that Eq.(32) 
reduces to Eq.(18) in the limit of Ms — t- with the relation '^j^ii—i^j)^^^^^^ = hj:)^^, as 
it should be. The fermionic contribution ^fermion " W ^f) can be obtained by replacing 
the factor -2 in V^^^^ ""^(a; Mg) by +21^/2] _ follows from the expression (32) that the 
contribution from a particle with the bulk mass Ms will be suppressed exponentially if 

Ms \. (37) 

This observation imply that Lifshitz particles of z > 1 even with Ms ^^/L may contribute 
to the effective potential if l/L y>> 1. This should be compared with the case of 2; = 1. 
Only massive particles with Ms ^^/L contribute to the effective potential. 

4 An SU{2) Gauge Model with an Adjoint Fermion 

Let us consider an SU{2) gauge theory on x Si coupled to a massless adjoint fermion 
whose kinetic term has the Hof ava type higher derivative. We are interested in the gauge 
symmetry breaking patterns and the mass of the adjoint scalar field which is originally 
the component gauge field for the compactified direction. 

The one-loop effective potential is useful tool in order to study them. In the present 
case, from the discussions in the previous section, the effective potential is given by 



^^^'(q;;M = 0) 



ir(f)^ - 



712 



3 £ 4 (1 + cos(27rn(2a))') 4hD=5,z 



n=l 



+ ( T ) ^wE^ > . ZilTT { 1 + cos(27rrz(2a)) 



L J r(5) ^ n 

' ^ ' n=l 



(38) 



The first hue in Eq. (38) is the contribution from the gauge (and ghost) fields and the 
second one is the one from the adjoint fermion. The Wilson hue phase a is related by the 
vacuum expectation value (Ay), 

gL{Ay) = 27r diag.(Q;, -a) (39) 

and the factor hD=5,z is given, from (19), by 

hD=5,z = z. (40) 
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The L is the length of the circumference of the circle S^. We see that the magnitude 
of the higher derivative, which is shown by the ratio l/L, plays the role of changing the 
size of the contribution from the adjoint fermion to the effective potential. The gauge 
symmetry breaking patterns through the Wilson line phases (Hosotani mechanism) has 
been studied extensively [27], and we understand that the pattern depends on the matter 
content of the theory. 

If the scale / is set to zero, then, the one-loop effective potential has the contribution 
from the gauge (and the ghost) fields alone. It has been known that the minimum of the 
potential is located at 

a = 1 mod - 
V 2 

for which the SU (2) gauge symmetry is not broken. On the other hand, if we take the scale 
/ to be large enough, then, the effective potential is dominated by the adjoint fermion. In 
this case, the minimum of the effective potential is 

a = 0.25, (42) 

so that the SU{2) gauge symmetry is broken to U{1). 

The above observation suggests that the gauge symmetry breaking patterns depend 
on the ratio l/L. Actually, the vacuum structure changes according to the magnitude of 



Veff 




Figure 1: The behavior of the one- loop effective potential with respect to the ratio = 
0.12,0.122,0.12237,0.123 with z = 2. According to the change of the ratio, the vacuum 
configuration changes from a = (mod 1/2) to a = 0.25. The degenerate vacuum appears 
at the critical ratio (-^)c = 0.12237, above which the vacuum configuration is a = 0.25, 
where gauge symmetry is broken down to U{1). 




l/L, as shown in Fig.l. We find that for z = 2 the degenerate vacuum appear at the 
critical ratio, 

= 0.12237. (43) 

For larger value than the critical ratio, the SU{2) gauge symmetry is broken to U{1), 
while the smaller value than the critical value, the SU{2) gauge symmetry is unbroken. 
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The behavior of the vacuum expectation value (VEV) a with respect to IjL is depicted 
in Fig.2, where we observe that the VEV jumps at the critical ratio, hence the phase 
transition is first order. 

a 
0.5- 

0.4- 

0.3- 

0.2- 

0.1- 

0.00 0.05 0.10 0.15 0.20 0.25^^^ 

Figure 2: The behavior of the vacuum expectation value (VEV) a with respect to At 
the critical ratio = 0.12237, the VEV jumps, and the phase transition is first order. 



Let us next discuss the mass of the adjoint scalar field which is originally the component 
gauge field for the direction. The mass of the adjoint Higgs scalar is obtained by the 
second derivative of the effective potential evaluated at the vacuum, 



2 r.OTrR'^XSl , X 



3^ 



7,2 



27r2L2 



oo ^ 

3 — cos(27m(2a)) 



>- n=l 



-4z 



r(4z + 1) f I 



L 



4(z-l) oo 



X:^cos(2vrn(2«)) 



n=l 



(44) 



where we have defined the 4d gauge coupling constant g = g/\/L. Hereafter we take 
z = 2 for numerical studies. 

Let us study the case of a = 0.25, for which the gauge symmetry is broken to U{1). 
The adjoint Higgs mass is given by 



9 , , 8! 
--C(3)+4x2x- 



27r2L2 

where we have used the Riemann's zeta function 



E 

n=l 



:-ii 



n 



4k-l 



-(1-2 



2-4fc 



)C(4A: - 1) 



(45) 



(46) 



^Let us note that the adjoint Higgs scalar can be massive even if the gauge symmetry is not broken. 
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with k — 1,2. The mass squared mjj is positive definite for the range of {l/L)c < l/L, 
where the vacuum configuration is given by a = 0.25. We note that the Gamma function 
r{4:Z + 1) is sizable enhancement for the mass. The adjoint scalar mass depends on the 
ratio l/L, which can also enhance the mass of the adjoint scalar field. 

The ratio between the adjoint Higgs scalar and the lightest massive gauge boson is 
given by 



rriH 



(0) 

w 
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4C(3) + 13230 X C(7) X ( - 




(47) 



where we have used the expression for the the massive Kaluza-Klein gauge bosons after 
the gauge symmetry breaking, 



(n)2 



W 



2n 

T 



{n + 2ay 

a=0.25 

We obtain the mass of the adjoint Higgs scalar field for various values of l/L 



for n = 0,1,2, 



(48) 



rriH 



(0) 



w 



' 0.54 


for 


/ 

L 


= 0.2, 


2.28 


for 


I 

L 


= 0.4, 


< 9.17 


for 


I 

L 


= 0.8, 


14.3 


for 


I 

L 


= 1.0, 


20.6 


for 


I 

L 


= 1.2, 



(49) 



The adjoint Higgs mass is generated through loop effects and is natural to be light com- 
pared with the massive gauge boson appeared after the gauge symmetry breaking. In the 
present case, however, thanks to the arbitrary scale I and the Gamma function V{A:Z-\- 1), 
the adjoint scalar mass can be heavier than the massive gauge boson, as shown above. 

The Higgs mass is also affected by the dynamical critical exponent z. In order to see 
it, let us study the ratio between the adjoint Higgs scalar and the lightest massive gauge 
boson with respect to z for fixed value of l/L. We take l/L — 0.25 as an illustration, for 
which the SU{2) gauge symmetry is broken to U{1). We obtain, from (44) and (48), that 



rriH 



(0) 



^ X <^ 



w 



0.12 


for 


z 


= 1, 


0.87 


for 


z 


= 2, 


7.50 


for 


z 


= 3, 


113.15 


for 


z 


= 4. 



(50) 



We observe that the Higgs mass highly depends on z and is enhanced by the effect of the 
higher derivative, as pointed out in the previous section. 

^If we take l/L ~ 0.393 (0.4) and g ~ 0.65, we have mu ~ 115 (119) GeV. 
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5 Higher Dimensional Extension 



In this section, we extend the previous analysis in section 3 to the Lifshitz type gauge 
theory on R^~^ x Sl^ x Sl^^ and show that the effective potential on R^'"^ x Sl^ x S].^ 
is given by the sum of an effective potential on R^~^ x R\^ x S^^ and infinitely many 
effective potentials for Kaluza-Klein modes on R^~^ x Sl^. We here focus on a massless 
scalar contribution to the one-loop effective potential 



scalar 



(ai,a2; Ms = 0) 



exp< —t 



mi=— oo 



m2=— oo 



X tr 



^2(21-1) ^ 27r(mi + Qi) j ^ ^2(22-1) /'27r(m2 + ^2) 



, (51) 



where the scalar action has been assumed to include ll^^^~^\Dy^y^^ and ll^^^~^\Dy^)^^^ 
for 5"^^ and Sl^ directions, respectively. The Li and L2 are the circumferences of Sl^ and 
Sl^, and Oil and a2 are defined by cti = g'(74j^^)Li/(27r) and a2 — g{Ay^)L2/{27r). We note 
that ai and 0:2 are assumed to commute each other to minimize the tree level potential 
5'^tr(adj)[^2/i, ^2/2]^- Using the Poisson summation formula (14), we obtain the UV finite 
expression 



scalar 



1 

4:71 



X tr(7^) 



(ai,a2;M5 = 0) 

n — 9 ^ 00 



L1L2 JQ 

m,Ti2=— 00 



/oo /"OO 
di^i / du 
-00 J —00 

exp{ -t [/^ (^) (^) 



— 27ri(niiii + 7121*2) + 27ri(niai + n2a2) | 



(52) 



where the prime of the summation denotes that the contribution of ni = ^2 = has to 
be removed from the summation to subtract the UV divergent part from the effective 
potential. The summation over ni and 712 may be rearranged as follows: 

00 00 00 

= (l-(^niA2,o) 
i=— 00 n2=— 00 
00 00 

^ Y Y {'^ni,o(l - (^n2,o) + (1 - (^ni,o)} 

ni=— 00 n2=— 00 

00 00 00 00 

= E E'+ E' E ■ (53) 



m,Ti2=— 00 ni=— oon2=— 00 
00 00 



n2=— 00 ni=— 00 n2=— 00 
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Thus, we can rewrite Eq.(52) as 

K^-^xsi xs: 



scalar 



(ai,a2; M = 0) 



n2=— oo 



X tr 



exp< 



1 

47r 



2 i 



LIL2 



ni=— 00 712=— 00 



X tr(7^) 



- 2TTi(niUi + ^2^2) + 2'Ki(niai + n2a2) | 



(54) 



The first term on the r.h.s. of Eq.(54) can be expressed, with the change of variable 
Pi — 2Trui/Li, as 



712 =— 00 



V47r/ 



dpi 
27r 



X tr(7e) 



duo 



27rin2U2 + 2'Kin20L2 



(55) 



It turns out that V^^aiar ('^2; M = 0) has a clear geometrical meaning. It is nothing 

but the one-loop effective potential of a D-dimensional Lifshitz type gauge theory on 
R^~'^ X X Sl^ (but not R^~'^ x 5*]^ x S]^) where a Lifshitz type higher derivative 
labeled by zi appears for one of the coordinates in R^~^ (i.e. R\^) as well as for the 
coordinate of Sl^ . 

The second term on the r.h.s. of Eq.(54) can be expressed, by using the Poisson 
summation formula reversely, as 



^ 00 



-1 



.2 — (47r)(^-2)/2L 

m2=— cx) ^ ' ni=— 00 



1 °o „oo 

- y^' dtt-''/^ / diii 



X tr(7e) 



scalar 



(56) 



m2=— 00 
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where 



M — 72:2-1 

^''-l-m.n <" 



^2 



27r(m2 + aa) ^2 ^^^^ 



Again, V^p^iar ('^i! ^m2) turns out to have a clear geometrical meaning. It is nothing 
but the one-loop effective potential of a (D — l)-dimensional Lifshitz type gauge theory 
(but not D-dimensional one) on x 5*^^ where the scalar action contains the Lifshitz 

type higher derivative of Zi for the 5"^^ direction with the bulk mass that corresponds 
to the Kaluza-Klein mass of the mode ma- Thus, we found that^'^ 

kL/ " "(«i=«2;M = 0) 

= Krr'''^^''^'^(«2;M = 0) + -l (58) 

^ m2=— 00 

Since the above observation is expected to hold for fermion and gauge fields, the one-loop 
effective potential of the Lifshitz type gauge theory on R^'"^ x 5"]^ x 5'^^ turns out be 
written into a decomposition form similar to Eq.(58). We also expect that the effective 
potential for a Lifshitz type gauge theory on R^~^ xSl^x- ■ - x S]^ has a similar structure, 
though we will not proceed furthermore. 



6 Conclusions and Discussions 

We have investigated the Lifshitz type gauge theory in the gauge- Higgs unification. The 
Lifshitz scalar and fermion possess the kinetic terms of the higher derivatives labeled 
by z for the direction of the extra dimension. We have succeeded to evaluate the one- 
loop effective potential for the zero mode of the extra dimensional component of the 
gauge field and found that it heavily depends on the dynamical critical exponent z of the 
Lifshitz particles. The overall sign of the effective potential can change with respect to 
z irrespective of the spin and statistics, and furthermore the one-loop effective potential 
turns out to vanish for even D and z. A huge numerical factor arises for 2; > 1 and this 
property may solve the light Higgs mass problem. It has also been found that the length 
parameter I plays an important role in determining the vacuum expectation value of the 
Higgs field as well as the magnitude of the Higgs mass. 

To obtain the finite expression of the one-loop effective potential, we have extended 
the space-time dimension D to a real (or complex) number and taken it to be the original 
integral value only at the final stage. To confirm this procedure, we have provided the 
mathematical tools and derived the one-loop effective potential in the several different 

similar structure has been found in gauge theories on R^~^ x at finite temperature [11]. 
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ways. We have also studied the one-loop effective potentials numerically in the 5d SU (2) 
model, and confirmed various peculiar properties of our models. 

The analysis for the Lifshitz type gauge theory on R^"^ x Si has been extended to 
the higher extra dimensions Sl^ x Sj^. An interesting observation is that the one- loop 
effective potential on R^""^ x Sj^ x Sl^ can be expressed by the sum of a one-loop effective 
potential on x Rl^ x Sl^ and infinitely many one-loop ones on R^~'^ x Sl^, which 

is one-dimension lower than the original space-time dimension D, coming from Lifshitz 
Kaluza-Klein particles on Sl^ . A similar decomposition property will hold for the Lifshitz 
type gauge theory on the higher extra dimension Sl^ x Sl^ x • ■ ■ x Sl^. 

Another extension is to replace the circle Si with the orbifold SI/Z2. We can apply 
the methods developed here to this case. Hence, it is important and interesting to study 
the gauge symmetry breaking and Higgs mass in a context of Lifshitz type gauge theories 
on SI/Z2. 

In this paper, we have assumed the higher derivative terms to present only for the 
direction of the extra dimensions. Hofava's original idea [18] is, however, to demand the 
anisotropy between time and space coordinates to make gravity theory power-counting 
renormalizable. According to the Hofava's spirit, we may replace the second order deriva- 
tives of all the spatial coordinates of R^~^ x S^ by higher derivatives to make the theory 
renormalizable. The power-counting renormalizability requires the mass dimension of the 
gauge coupling to be non-negative. This implies that the dynamical critical exponent z 
should be greater than or equal to D — 3. The one- loop effective potential coming from 
such a Hofava- Lifshitz massless scalar loop may be given by 



where 



scalar 



a;M = 0) 

^ ' m=— 00 " * 

r ,9r. n// o\z /27r(m + q;)\ 2^ 



X tr 



sm{n{z + D- 2)/2)r(^ + D - 2) 



n=l 



LD\L 



sui{ti{z + D- 2)/{2z))r{{z + D- 2)/{2z)) 
cos(27ma) 



n 



z+D-l 



(59) 



(60) 



Wc should emphasize that standard gauge-Higgs unification models are not renormalizable 
because the total space-time dimension D is greater than 4. Thus, only a very limited 
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class of physical quantities such as the Higgs mass are finite and calculable, so that the 
predictability of the theory is quite restricted. Since the Hof ava-Lifshitz type gauge theory 
is power-counting renormalizable, any physical quantities can, in principle, be computed 
with finite values at the cost of Lorentz symmetry violation. It would be of great interest 
to investigate Hof ava-Lifshitz type gauge theory in more details. 
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A Other Methods to Evaluate the Effective Potential 

In the foUowings, we show the alternative ways to reproduce the result (18). 
A.l A method with contour Integrals 

Eq. (18) can be reproduced in a way similar to [38]. First, in terms of contour integrals, 
Eq. (22) can be written as 

KSa"'^^'(«;^5 = o) 



where 



(61) 



F{w,a) — cos{2mv) — cos{2'Ka), (62) 



and the contour C is the set of the circles each of which surrounds w = m ± a > 0, 
m = 0, 1, ■ ■ • , on the plane. Since the integrand is regular except for roots on the real 
axes, we can change a closed path of Re w > to the contour shown in Fig. 3, where the 
closed path consists of the infinite semicircle on the plane Ci, the straight lines C2 from 



18 




► Re(w) 



Figure 3: The integration path stretched on the positive-real plane. 

-|-ioo to +i5 {5 is a positive infinitesimal) and C4 from —i5 to —ioo, and a semicircle C3 
centered at the origin with radius 6 . The integrals along the infinite semicircle vanishes 
in the sense of analytic continuation and by taking the limit (5^0, only integrals along 
C2 and C4 remain non-vanishing, and thus we obtain 



scala. («;M5 = 0) 

r (-^) ^27rV^''"'^ z{D - 1) 

L 



(47r) 2 



27ri 



xtr 



(^) 
+ 



/■oo 

1) / (iU^(^-^)-Mn[F(-2t, a)] 
^0 



r (-^) /27r\ ^(^-'^ - 1) 



(47r)'^ ^ 
'7rz(L' - 1; 



X sm 



L 



TT 



where we have used F{w, a) — F{—w, a) and 

poo 

Jo 



Recalling a relation cosh(A) — cos(S) = e^{l — e 
as 



-A+iB 



)(1 



(63) 
(64) 

"^~*^)/2, we can rewrite J 



J ^ dt {ln(l - e-2'^*+^2,ra) + in(i _ ^-2^t-i2^a^ ^ 27ri - ln(2) } . (65) 

Jo 
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We note that the above integral contains divergences. We ignore such divergent terms 
because we are interested in a-dependent ones. Using the expansion formula 



ln(l - e-'^*+^2^"") = ' (66) 

n=l ^ 

and after integration we get 

_ 2 , , cos(27rnQ;) , , 

J = (67) 

^ ^ n=l 

Combining (63) and (67), and utilizing the Gamma function formulas, one can reproduce 
the same result as (18). 

A. 2 Use of Hypergeometric Functions 

Here we show the one more way to evaluate the effective potential. Using Poisson sum- 
mation, we can write (8) as 



v^(27r)^(-°-i) 1 fl'-^^""-' 



D-l 



oo 



tr(^) J2 e'''^"" / dr T-^-'F,, (27rn; r) , (68) 



where 

1 f°° 

F2,{y;T) = -= / dx exp[-Tx'^'yy (69) 
v27r J-oo 

and we have used the following Poisson summation formula 

oo oc 

^ exp[-T(m + «)'"] = V2^ ^ e^^'^^'^Fs, (27rn; r) . (70) 

m=— oo n=— oo 

We can rewrite F2z{y; r) into the form 

\ 2r 

ly 



1 ^ I rV(2^) J /2r + l\ 
^"'''•"'° V2^.rV.^^) S (2r)! ) ' 

where we have expanded e'^^ in powers of x and used the integral representation of the 
gamma function. To proceed further, we note that any non-negative integer r can uniquely 
be parameterized by two integers s and k such that r = zs + k, where r = 0,1,2,---, 
s = 0, 1, 2, • • • and A; = 0, 1, • • • ,z — 1. In terms of the Pochhammer's symbol 

(a), = a{a + l){a + 2)...{a + s-l)^ ^^^^ , (72) 
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we have 



1 ^ r(g|±l) / ty y 



It then foUows that^^ 

A;=0 

where we have used the relation (l)^ = s\ and pFq[- • • ) is the generahzed hypergeometric 
function defined by 

■^-^ {0.1) s{0-2) s ■ ■ ■ {o-p) s 

pFg{ai, • • • , Op; 61, • • • , 6,; w) = 2^ (h ^ (h ^ TT^TT • i^^) 

Truncating the n — mode, which is independent to a and responsible for the UV 
divergence, and integrating with respect to r, we obtain 



2\/27rf27rl^(^-i) 1 /p-^X^"^ ^ 
= - ; L 7 W tr(^)E^°^(2^^«)^^^(27rrz,i^), (77) 
(47r) 2 L \ L J 



(47r) 

where G2z{y, D) is defined by 

/■oo 

G2.(i/,i^) = / dTT-'^--^F2,{y;T). (78) 

By use of the formulas of an indefinite integral and an asymptotic expansion^^ of gener- 
alized hypergeometric functions 

J du m""-^ pFq{ai, ■■ ■ , Op; 61, • • • , bg] u) 

= — p+iFq+i{a,ai, ■ ■■ ,ap;a + l,bi, - ■ ■ ,bq;u), (79) 
pFq{ai, ■ ■ ■ , Op; 61, ■ • • , bq] u) 

,„j_ ^ rK)(n,t.rfe))n;..,,,r(a,-«,) ^ 

(n?.,r(«,))nj.,m-ao 



^^The authors would like to thank Professor H. So who taught us the derivation of Eq. (75) given here. 

"^^For z > 2, in the asymptotic expansions we have extra terms, which are exponentially diverging and 
may spoil the finiteness of G2ziy,D) (and therefore the a-dependent part of the effective potential may 
suffer from divergences). However, by numerical studies we found that such divergences in G2z{y,D) 
are canceled out with each other in different fc's and finite values of Gzziy-D) are obtained for smaller 
z's. Therefore we have omitted such terms in (80). See Ref. [39] for detailed asymptotic expansions of 
generalized hypergeometric functions. 
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we can show after some calculations that 

2D-1/2 Y{D/2)T{z{D - 1) + 1) sm(^^7r) 



y 



z{D-l)+l 



r{D) 



CD-I), 



sin(^-^7r) 



Inserting it into Eq. (77) leads to Eq. (18), as expected. 



B Effective Potential for Massive Matter 



In this appendix we show the calculation for (32) in detail. Let us start with 

^{a, M) 

\ / m — — rv^ 



2(0F) 



D-l 



(27r)^-ir(^) L 



X 



J2 / dEE{VE'' - M'^f-hr^n) In 



a I ]2z-2 



27r(m + a) 
L 



2z 



For odd D, we can write 



where „Cfe = n!/(A;!(n — A;)!) is the binomial coefficient. Iq is given by 



L = 



roc 



dEE'^'i+hiiji) In 



m=— 00 

_1 N 2g+2 „oo 



^2 ^ ;2z-2 



27r(m + «) 



2z 



ML'=/l^-'^ 



where 



W= ln[g' + (27r(m + a))2^]. 

m=—oo 

Here we have ignored a- independent terms. Since the derivative dW/ dQ is given by 



dW 
'dQ 



E 



2Q 



G',{Q) 



2z 



G,(Q) ^ sin (^a - ^ij|-^) 



u!j = exp 



^{2J - l)7r 
2^ 
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(the derivation of (86) is given in the latter part of this appendix), we obtain 

2z 



W ^ In 



In 



= In 



■ z 

u 

■j=l 
1 

4 



sm TT a — 



27r 



^ ^ sin ^TT — 



27r 



1 ^ 

- J] exp(s,Qi/^) {1 - exp[-Qi/^s,- + i{27ra - Q'^'cj)]} 



X 



{1 - exp[-g^/^Sj- - i{27Ta - Q^^'cj)]} 



(89) 



where we have used the relation ouj — uj2z+i-j, and 



Imujj — sin 



Reujj — cos 



2z 

^(2J• - l)7r 
2z 



(90) 
(91) 



Neglecting some a-independent parts in W which correspond to the ultraviolet diver- 
gences, we write Ig as 

_1 \ 2q+2 z 



\ / j=l 



poo 

Ja,j = / dQQ^'^+^]n[l-exp{-{sj + icj)Q^/' + i2Tia)]. 



Using (66) and performing the integration of Jgj, we get 



(92) 
(93) 



zySj + iCj^ 



'z{2q+2) 



E;^i(^^r(.(2, + 2),np,) 



(94) 



n=l 



where T{x,y) = ^dx is the incomplete Gamma function and 



Combining (83), (92) and (94), we obtain 



scalar 



a, M) 



D-3 

2 



2(0F)^-^ 1 >p ^, ^2. 



^-q I I' 

2 



1 \ 2g+2 



°° gi27rna 



2:(2q+2)+l 



n=l 



n 



r{z{2q + 2), 



2z 

in LM^/' 



(95) 
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One can rewrite (95) into more useful form. Utilizing an expansion formula r(n, x) — 
{n — l)!e~^ x'^/A;!, we can rewrite Jqj as 



n=l k=0 
zi2q+2)-l 

= -^[^(2g + 2)-l]!(-za;,)^(^^+^) J] ^ Li,(,,+,)_,+i(e-''^+^2™), 



k=0 



(96) 



where we have changed the order of the summation and Lis{x) is defined in (35). Com- 
bining (83), (92) and (96), we obtain (32). 

Here we outline the derivation of (86). At first we point out that the left-hand-side of 
(86) can be written in terms of the following contour integral: 



E 



(27r(m-Fa))2^ 



dw 



TTCOtfTTw) 



Q (27r(w + a))2^' 



(97) 



where the path Ci denotes the set of the circles surrounding the points w = m. Without 
changing the value of the integral, we can replace the integration path Ci by Cii+ -|- Cii_ 
shown in Fig. 4. Counting the residues of the integrand at the points which are enclosed 




► Re(w) 



Figure 4: Two contours Cii±. Cii+ (Cii-) consists of the infinite semicircle above (below) 
the line Im(w) = 0, and the straight path from -|-oo -|-ie (— oo — ie) to —oo + ie {+oo — ie). 
e denotes a positive infinitesimal. 
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by Cii±, we obtain 

TT cot(7rw) 

= -t , (98, 

fc=l,fc^j 

where flj = {xcuj — 27rQ;)/27r, and ReStt,_>a[/(tt;)] denotes the residue of f{w) at w — a. 
One can reproduce (86) by plugging x = Q^^^ to (98) and using the formula 

n (uj-Uk)^--. (99) 

References 

[1] M. S. Manton, Nucl. Phys. B158, 141 (1979). 

[2] D. B. Fairlie, Phys. Lett. B82, 97 (1979). 

[3] J. Sherk and J. Shwarz, Phys. Lett. B82, 60 (1979); Nucl. Phys. B153, 61 (1979). 

[4] Y. Hosotani, Phys. Lett. B126, 309 (1983), Ann. Phys. (N.Y.) 190, 233 (1989). 

[5] C. Csaki, C. Grojean, H. Murayama, L. Pilo and J. Terning, Phys. Rev. D69, 055006 
(2004). 

[6] M. Sakamoto, M. Tachibana and K. Takenaga, Phys. Lett. B458, 231 (1999); Prog. 
Theor. Phys. 104, 633 (2000). 

[7] M. Sakamoto, M. Tachibana and K. Takenaga, Phys. Lett. B457, 33 (1999). 

[8] S. Matsumoto, M. Sakamoto and S. Tanimura, Phys. Lett. B518, 163 (2001); M. 
Sakamoto and S. Tanimura, Phys. Rev. D65, 065004 (2004). 

[9] H. Hatanaka, K. Ohnishi,M. Sakamoto and K. Takenaga, Prog. Theor. Phys. 107, 
1191 (2002), Prog. Theor. Phys. 110, 791 (2003). 

[10] K. Ohnishi and M. Sakamoto, Phys. Lett. B486, 179 (2000); H. Hatanaka, S. Mat- 
sumoto, K. Ohnishi and M. Sakamoto, Phys. Rev. D63, 105003 (2001). 

[11] M. Sakamoto and K. Takenaga, Phys. Rev. D80, 085016 (2009). 



- — - Q) dw 



vrcotfTTw) 



27ri Jc, x^^ + (27r(w + a)f^ 



= -E, 



Res 



25 



[12] L. Randall and R. Sundrum, Phys. Rev. Lett 83, 4690 (1999). 

[13] H. Hatanaka, M. Sakamoto, M. Tachibana and K. Takenaga, Prog. Theor. Phys. 102, 
1213 (1999). 

[14] T. Nagasawa and M. Sakamoto, Prog. Theor. Phys. 112, 629 (2004). 

[15] M. Sakamoto and K. Takenaga, Phys. Rev. D75, 045015 (2007). 

[16] L. Hall, Y. Nomura and D. Smith, Nucl. Phys. B639, 307 (2002); Y. Hosotani, 
S. Noda and K. Takenaga, Phys. Lett. B607, 276 (2005). 

[17] K. Agashe, R. Contino and A. Pomarol, Nucl. Phys. B 719, 165 (2005); Y. Hosotani 
and M. Mabe, Phys. Lett. B615, 257 (2005); K. Agashe and R. Contino, Nucl. Phys. 
B 742, 59 (2006); K. y. Oda and A. Weiler, Phys. Lett. B 606, 408 (2005); Phys. 
Rev. D 73, 096006 (2006); R. Contino, L. Da Rold and A. Pomarol, Phys. Rev. D 
75, 055014 (2007); A. Falkowski, Phys. Rev. D 75, 025017 (2007); A. D. Medina, 
N. R. Shah and C. E. M. Wagner, Phys. Rev. D 76, 095010 (2007); H. Hatanaka, 
arXiv:0712.1334 [hep-th]. N. Haba, S. Matsumoto, N. Okada and T. Yamashita, Prog. 
Theor. Phys. 120, 77 (2008); C. Csaki, A. Falkowski and A. Weiler, JHEP 0809, 008 
(2008). 

[18] P. Hofava, Phys. Rev. D79, 084008 (2009). 

[19] T. Takahashi and J. Soda, Phys. Rev. Lett. 102, 231301 (2009); G. Calcagni, JHEP 
0909, 112 (2009); E. Kiritsis and G. Kofinas, Nucl. Phys. B821, 467 (2009); S. 
Mukohyama, J CAP, 0906, 001 (2009); R. Brandenberger, Phys. Rev. D80, 043516 
(1909); Y.-S. Piao, Phys. Lett. B681, 1 (2009); S. Mukohyama, K. Nakayama, F. 
Takahashi and S. Yokoyama, Phys. Lett. B679, 6 (2009); S.K. Rama, Phys. Rev. 
D79, 124031 (2009); R.A. Konoplya, Phys. Lett. B679, 499 (2009); B. Chen, S. 
Pi and J.-Z. Tang, JCAP, 0908, 007 (2009); E.N. Saridakis, Eur.Phys.J. C67, 229 
(2010). 

[20] R.-G. Cai, L.-M. Cao and N. Ohta, Phys. Rev. D80, 024003 (2009); R.-G. Gai, Y. Liu 
and Y.-W. Sun, JHEP, 0906, 010 (2009); Y.S. Myung and Y.-W. Kim, Eur.Phys.J. 
C68, 265 (2010); A. Kehagias and K. Sfetsos, Phys. Lett. B678, 123 (2009); R.-G. 
Gai, L.-M. Gao and N. Ohta, Phys. Lett. B679, 504 (2009); Y.S. Myung, Phys. Lett. 
B678, 127 (2009); R.B. Mann, JHEP, 0906, 075 (2009); S. Chen and J. Jing, Phys. 
Lett. B687, 124 (2010); S. Ghen and J. Jing, Phys. Rev. D80, 024036 (2009). 

[21] H. Lu, J. Mci and C.N. Pope, Phys. Rev. Lett. 103, 091301 (2009); H. Nikolic, Mod. 
Phys. Lett. A25 1595 (2010); E.G. Colgain and H. Yavartanoo, JHEP 0908 (20021) 



26 



09; T. Sotiriou, M. Visser and S. Weinfurtner, Phys. Rev. Lett. 102, 251601 (2009); 
R.-G. Cai, B. Hu and H.-B Zhang, Phys. Rev. D80, 041501 (2009); D. Orlando and 
S. Reffert, Class. Quant. Gmv. 26, 155021 (2009); T. Nishioka, Class. Quant. Gmv. 
26, 242001 (2009); C. Charmousis, G. Niz, A. Padilla and RM. Saffin, JHEP 0908 
(20070) 09; M. Li and Y. Pang, JHEP 0908 (20015) 09; J. Chen and Y. Wang, 
Int. J. Mod. Phys. A25, 1439 (2010); T.R Sotiriou, M. Visser and S. Weinfurtner, 
JHEP 0910 (20033) 09; Y.-W. Kim, H.W. Lee and Y.S. Myung, Phys. Lett. B682, 
246 (2009); M. Sakamoto, Phys. Rev. D79, 124038 (2009). 

[22] D. Anselmi and M. Halat, Phys. Rev. D76, 125011 (2007); R. Icngo, J.G. Russo and 
M. Serene, JHEP 0911 (20020) 09; M. Visser, Phys. Rev. D80, 025011 (2009); D. 
Anselmi, Ann. Phys. 324, 874 (2009); D. Anselmi, Ann. Phys. 324, 1058 (2009); D. 
Orlando and S. Reffert, Phys. Lett. B683, 62 (2010). 

[23] S.R. Das and G. Murthy, Phys. Rev. D80, 065006 (2009); A. Dhar, G. Mandal and R. 
Wadia, Phys. Rev. D80, 105018 (2009); A. Dhar, G. Mandal and P. Nag, Phys. Rev. 
D81, 085005 (2010); J, Alexandre, F. Farakos, P. Pasipoularides and A. Tsapahs, 
Phys. Rev. D81, 045002 (2010); J, Anagnostopoulos, F. Farakos, P. Pasipoularides 
and A. Tsapalis, arXiv: 1007.0355 [hep-th]; 

[24] S.R. Das and G. Murthy, Phys. Rev. Lett. 104, 181601 (2010); M. Gomes, et al, 
Phys. Rev. D81, 045013 (2010); C.M. Reyes, Phys. Rev. D82, 125036 (2010); J.M. 
Romero, J. A. Santiago, O. Gonzalcz-Gaxiola and A. Zamora, Mod. Phys. Lett. A25 
3381 (2010); J. A. Alexandre, arXiv: 1009.5834 [hep-ph]. 

[25] A.F. Ferrari, H.O. Girotti, M. Gomes, A.Yu. Petrov and A.J.da Silva, 
arXiv:1006.1635 [hep-th]. 

[26] G. Panico, M. Serone and A. Wulzer, Nucl. Phys. B739, 186-207 (2006). 

[27] A. T. Davies and A. McLachlan, Nucl. Phys. B317, 237 (1989), A. McLachlan, Nucl. 
Phys. B338, 188 (1990), J. E. Hetrick and C. L. Ho, Phys. Rev. D40, 4085 (1989), 
C. L. Ho and Y. Hosotani, Nucl. Phys. B345, 445 (1990), A. McLachlan, Nucl. Phys. 
B338, 188 (1990), H. Hatanaka, Prog. Theor. Phys. 102, 407 (1999), K. Takenaga, 
Phys. Lett. B425, 114 (1998); Phys. Rev. D58, 026004 (1998); 66 085009 (2002); 
Phys. Lett. B570, 244 (2003); N. Haba, K. Takenaga and T. Yamashita, Phys. Lett. 
B605, 355 (2005). 

[28] B. Chen and Q-G. Huang, Phys. Lett. B683, 108 (2010); P. Hofava, Phys. Lett. 
B694, 172 (2010). 



27 



[29] M. Gomes, J-R. Nascimento, A.Yu. Petrov and A.J.da Silva, Phys. Lett B682, 229 
(2009). 

[30] W. Xue, arXiv:1008.5102 [hep-th]. 

[31] Y. Kawamura, Prog. Theor. Phys. 122, 831 (2009); K. Kaneta and Y. Kawamura, 
arXiv:0909.2920 [hep-ph]; W. Chao, arXiv:0911.4709 [hep-th]; D. Anselmi, Phys. 
Rev. D79, 025017 (2009); D. Anselmi, Eur. Phys. J. C65, 523 (2010); M. Pospelov 
and Y. Shang, arXiv: 1010.5249 [hep-th]; D. Anselmi and E. Ciuffoh, Phys. Rev. D83, 
056005 (2011). 

[32] A. Padilla, J.Phys.Conf.Ser. 259, 012033 (2010); D.Blas, O. Pujolas and S. 
Sibiryakov, JHEP 1104 (2011) 018; T. P. Sotirou, J.Phys.Conf.Ser. 283, 012034 
(2011). 

[33] E. M. Lifshitz, Zh. Eksp. Theo. Fiz. 11, 255 (1941), ibid. 11, 269 (1941). 

[34] M. Chaikin and T. C. Lubcnsky, Principles of Condensed Matter Physics, (Cam- 
bridge University Press, Cambridge, 2000). 

[35] S. Moriguchi, K. Udagawa and S. Hitotsumatsu, Iwanami Sugaku Koushiki 11, 
(Iwanami Press, Tokyo, 2000). 

[36] I.S. Gradshteyn and I.M. Ryzhik, TABLE OF INTEGRALS, SERIES,, AND PROD- 
UCTS, (Academic Press, New York, 1980). 

[37] A. Delgado, A. Pomarol and M. Quiros, Phys. Rev. D60, 095008 (1999). 

[38] J. Garriga, O. Pujolas, T. Tanaka, Nucl. Phys. B605, 192-214 (2001). [hep- 
th/0004109]. 

[39] The Wolfram Functions Site (http://functions.wolfram.com/07.31.06.0041.01). 



28 



